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' Abstract 

O ' 

The effective string action of the color-electric flux tube in the U(l)xU(l) dual 
Ginzburg-Landau (DGL) theory is studied by performing a path-integral analy- 
Oh, sis by taking into account the finite thickness of the flux tube. The DGL theory, 

corresponding to the low-energy effective theory of Abelian-projected SU(3) gluo- 
dynamics, can be expressed as a [U(l)] 3 dual Abelian Higgs (DAH) model with a 
certain constraint in the Weyl symmetric formulation. This formulation allows us to 
adopt quite similar path-integral techniques as in the U(l) DAH model, and there- 
fore, the resulting effective string action in the U(l)xU(l) DGL theory has also 
quite a similar structure except the number of color degrees of freedom. A modified 
Yukawa interaction appears as a boundary contribution, which is completely due 
to the finite thickness of the flux tube, and is reduced into the ordinary Yukawa 
interaction in the deep type-II (London) Limit. 
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1 Introduction 



The dual superconducting picture of the QCD vacuum [1,2] provides an in- 
tuitive scenario on the quark confinement mechanism. The formation of the 
color-electric flux tube can be understood analogously to that of the mag- 
netic Abrikosov vortex [3,4,5] in an ordinary superconductor via the Meissner 
effect with the electric Cooper-pair condensation. In fact, in the QCD vac- 
uum, color-magnetic monopoles are condensed, and then, the color-electric 
field associated with quarks is squeezed into an almost one-dimensional ob- 
ject forming the color-electric flux tube due to the dual Meissner effect (the 
dual Higgs mechanism). The linear rising potential naturally appears between 
the quark and the antiquark, which explains why a quark cannot be isolated 
as a free particle. 

Remarkably, color-magnetic monopole condensation is observed by Monte- 
Carlo simulations of lattice gluodynamics in the so-called maximally Abelian 
(MA) gauge [6,7,8]. In the MA gauge, Abelian dominance [9,10] for the non- 
perturbative quantities, such as the string tension [11,12] and chiral conden- 
sate [13], is also observed. Moreover, it is suggested that monopole configura- 
tions lead to large contributions to the string tension [14] and chiral conden- 
sate [15]. These results indicate that Abelian components of the gluon field 
and monopoles become relevant degrees of freedom to describe the low-energy 
(long-range) properties of the vacuum. 

These numerical results have justified a posteriori a low-energy effective the- 
ory of the gluodynamics composed of Abelian degrees of freedom (Abelian- 
projected gluodynamics [16]) with monopole condensation. The resulting ef- 
fective action is known to have the form of the Ginzburg-Landau (GL) theory 
which in its original version would effectively describe bulk properties of the 
ordinary superconductor. Since now the role of electric and magnetic prop- 
erties is interchanged, this is called the dual Ginzburg-Landau (DGL) the- 
ory [17,18,19,20]. The DGL Lagrangian is composed of the dual gluon field 
interacting with the monopole field, and in fact, the color-electric flux-tube 
configuration is obtained as the classical solution when the test electric charges 
are put into the vacuum. The thickness of the flux tube is controlled by the 
inverse masses of the dual gauge boson and the monopole, that is, the penetra- 
tion depth nig 1 and the coherence length m' 1 . The latter length characterizes 
the so-called core size of the flux tube. 

Although long-distance aspects are dominated by Abelian-projected degrees 
of freedom, the DGL theory may not be a complete substitute for all aspects 
of low-energy gluodynamics, since there would be an interface in the real 
theory with short range physics which is taken care of by charged vector and 
ghost fields. Such effects would play an important role for the description of 
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a dynamical charge density around test charges. Nevertheless, it is interesting 
to remind that the Abelian-projected flux-tube profiles have also been studied 
within the lattice gluodynamics in the MA gauge [21,22,23,24,25,26], which 
have shown quite similar profiles as the classical flux-tube solution in the DGL 
theory. Therefore it will make sense to regard the DGL theory as a starting 
point of further investigations of the hadronic flux tube [27,28,29,30] with the 
pointlike approximation of color-electric test charges. P] 

From the Abelian-projected SU(2) gluodynamics one can construct the U(l) 
DGL theory (often called the U(l) dual Abelian Higgs (DAH) model), where 
the U(l) dual symmetry corresponds to the Cartan subgroup of SU(2). One 
can analyze the physical essence of the dual superconducting vacuum and, 
at the same time, the structure of the flux tube. For the description of the 
dynamics of the flux tube, it makes sense to use a framework which directly 
deals with the string-like properties of the flux tube. Such framework is the 
string representation or effective string action approach of the DGL theory, 
which can be obtained by performing a path-integral analysis. In the U(l) 
case, based for simplicity, on the London limit approximation m x — > oo , 
the string representation was derived in Refs. [31,32,33,34]. The London limit 
approximation corresponds to the neglect of the core size of the flux tube. The 
resulting leading action was shown to be the Nambu-Goto (NG) action with 
the rigidity term [35,36]. The string representation of the U(l)xU(l) DGL 
theory corresponding to the Abelian-projected SU(3) gluodynamics was also 
studied using this approximation [37,38]. 

Recently, we have considered the derivation of the string representation of the 
U(l) DGL theory beyond the London limit, i.e. we have taken into account 
the finite thickness of the flux-tube core by keeping the mass of the monopole 
finite [39]. The used path-integral analysis combined with differential form 
techniques turned out to provide a very compact and transparent method for 
studying these problems including, in particular, also boundary terms related 
to the non-confining part of the potential. As a main point, the finite core 
corrections to the effective string action have been derived. 

In this paper, we are going to generalize the above work by investigating the 
U(l)xU(l) DGL theory beyond the London limit. Note that two of us have 
shown that the Weyl symmetric formulation reduces the U(l)xU(l) DGL 
theory into the sum of three types of the U(l) DAH model, which was quite 
useful to simplify the treatment of the DGL theory and the investigation of 
the flux-tube solution [29,30]. Therefore, in order to derive the effective string 
action, we pay here special attention to the Weyl symmetric formulation when 
using again path-integral techniques in differential form notation. 



The resulting DGL description is thus, in particular, justified for the string away 
from the charges and their dynamical charge clouds. 
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This paper is organized as follows. In Sec. 2 we introduce the Weyl symmetric 
formulation of the DGL theory in differential form notation. In Sec. 3 the role 
of the open Dirac string with color-electric charge at its ends for the viola- 
tion of the dual Bianchi identity is considered. In Sec. 4 the functional path 
integration in differential form notation is performed which requires special 
attention to gauge conditions and the Weyl symmetry constraint. This leads 
to the effective string action, where the finite thickness of the core of the flux 
tube is taken into account. Finally, Sec. 5 contains the summary. 



2 Weyl symmetric formulation of the DGL theory 



The Weyl symmetric expression of the DGL action in the differential form is 
given by 



'DGL 



= E 

i=i 



^ + ((d 



l B t )xl ) {d + iB i ) Xi ) + \{\ Xt \ 2 -v 2 ) 2 



(1) 



where the ingredients of the DGL theory are three types of the dual gauge 
field Bi (J2i=iBi = 0, i = 1,2,3), and the complex-scalar monopole field 
Xi = 0iexp(i?7j) ((pi, rji G dt, i = 1,2,3). In the presence of external quark 
sources in the vacuum describing a static quark-antiquark (q-q) system, the 
dual field strength Fj is expressed as 



Fi = dBi - 2tt * £, 



(m) 



(2) 



with 



-,(m) 



= E 



m 



'■J 



(3) 



Here Tr™ (i = 1, 2, 3) denotes the "color-magnetic" (Weyl symmetric) repre- 
sentation of the color-electric Dirac string (j — 1, 2, 3) which describes the 
open color-electric Dirac string [30]. Let us briefly summarize the definitions 
of operators in differential forms and useful formulae in Table 1, and quote 
the ingredients of the DGL theory in Table 2. The factor 2iririij in front of the 
color-electric Dirac string appears as a result of the extended Dirac quantiza- 
tion condition between the color-magnetic charges = gti (i = 1, 2, 3) and 
color-electric charges = ewj (j = 1,2,3): 

Qi m) • Qf = 9?i ■ ewj = 2-Krriij {eg = 4n), (4) 
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Table 1 

Definitions in differential forms in four-dimensional Euclidean space-time. 



r-form (0 < r < 4) 


CO 


w = H^i-Mr^i A ... A dx^ r 


exterior derivative 


d 


r-form i— ► (r + l)-form 


Hodge star 


* 


r-form (4 — r)-form 




** 


multiply by a factor (— l) r for an r-form 


codifferential 


5 = — * d* 


r-form (r — l)-form 


Laplacian 


A = d5 + 5d 


r-form i— > r-form 


Inner product 


(lu,t]) = f lo A *rj 


(w, »?) = H / d 4 xuj lll ... flr r] fll ... flr (w, r? G r-form) 






(w) 2 = (w.w) 



Table 2 

Ingredients of the DGL theory in differential-form notation. The subscripts z = 
1, 2, 3 denote the label of color in the color- magnetic (Weyl symmetric) representa- 
tion, 

dual gauge field 1-form B, L = B{ ^dx^ 

monopole field 0-form Xi = 4>% exp(zr/i) 

electric Dirac string 2-form S-" 1 - 1 = ^E^dx^ A cfaj, 

electric current 1-form j[ m ^ = jf^ dx^ 

where q and Wj are the root and weight vectors of the SU(3) Lie algebra, 
respectively. Explicitly, these vectors are defined as t\ = (— 1/2, y/3/2), e*2 — 
(-1/2,-^/2), e 3 = (1,0) and m = (1/2,^/6), w 2 = (-1/2,^3/6), w 3 = 
(0,-1/^). The labels j = 1,2,3 describe, for instance, red (R), blue (B), 
and green (G) color-electric charges, respectively. The matrix 

3 

rriij = 2ei -Wj = Yl £ ijk (5) 
fc=i 

takes the integer values ±1 or 0. Based on these definitions of charges, the 
R-R system is defined by £j=i = E ^ 0, ^=2 = 0' anc ^ ^j=3 = 0> which is 
also expressed as = 0, = — an d £|=3 = E in the color-magnetic 
representation. Clearly, the color-electric Dirac string in the color-magnetic 
representation satisfies the relation Yh=i = ( see > Table 3). 

The physics described by the DGL theory is characterized by two mass scales, 
the masses of the dual gauge boson m# = ^2/ f3 g v = y/3gv and of the 

monopole m x = 2y r \v, where (3 g = 2/(3g 2 ) is the dual gauge coupling and A 
the self-coupling of the monopole field. Recall that their inverses correspond to 
the penetration depth and the coherence length, respectively. We can further 
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Table 3 

The color-electric Dirac string structures in the color-electric and color-magnetic 
representations for the q-q system, where E corresponds to the world sheet of one 
string singularity. 





color-electric 


rep. 


color- 


■magnetic 


rep. 






y(e) 




v (m) 


v (m) 


v (m) 


R-R 


s 











-E 


E 


B-B 





E 





E 





-S 


G-G 








E 


-E 


E 






define the Ginzburg-Landau (GL) parameter k as 



k = m x /m B , (6) 

which controls the vacuum property of dual superconductivity, in particular, 
whether the vacuum belongs to type-I (k < f) or type-II (k > 1) or to their 
border (k — 1) . 



3 Open string singularity and the violation of the dual Bianchi 
identity 

When we consider the open string system, such as the R-R system, it is 
important to pay attention to the color-electric Dirac string singularity in 
the dual field strength F i: since it is directly related to the violation of the 
dual Bianchi identity: 

+2 7 r* J fV0, (7) 
— * are used. Since the latter relation is 

(8) 

one finds that the color-electric current is defined as the boundary of the 
Dirac string T!f l \ In order to treat the Dirac string singularity further, it is 
useful to decompose the dual gauge field Bi into two parts called regular and 
singular ones, 

B t = BP + Br is , (9) 



dFi = d 2 Bi - 2nd * E\ m) 



where d 2 Bi = and d * E 
also written as 

xy M _ Am) 



(m) 
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where the singular part is explicitly written as 

J Bf ng = 27rA- 1 5*£ l (m) . (10) 

The inverse of the Laplacian, A -1 , is the Coulomb propagator (AA -1 = 
A -1 A = 1). Note that A and A -1 commute with the exterior derivative 
d, the co differential 5, and the Hodge star *. Notice further that _Bf ng is a 
background gauge field associated with the singular string surface. By using 
the relation dA5 + 5 Ad = 1, one finds that Eq. (10) leads to 

dBf ng = 27r(*sl m) + *Cf ) ). (11) 

Then, the dual field strength is written as 

F i = dBl es + 2ir*ct\ (12) 

where 

*ct ) = &- 1 8*jt ) - (13) 



Clearly, from this expression we see the violation of the dual Bianchi identity, 
since d 2 Bl° g = and d * C- m ^ = Then, the DGL action is written as 



3 

Sdgl = ^2 
i=i 



(14) 



where the polar decomposition form of the monopole field Xi = 0i ex P(^i) 
(<pi, r/i G 9ft) is inserted. Since we are interested in the oj^en-string system, 
we investigate the case that the phase of the monopole field is single- valued, 
so that we have i] = i] TCg . In general the phase r\ can contain the singular 
part, which describes a closed-string, corresponding physically to the glueball 
system [28]. Such contribution should be taken into account when glueballs 
and quantum effects in the dual superconducting vacuum are concerned. 



4 Functional path integration and the effective string action 

The effective string action S e g (S^), written in terms of the color-magnetic 
representation of the color-electric Dirac string, will be obtained after the 
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functional path-integration over both the dual gauge fields and the monopole 
fields. Explicitly, this can be expressed by the following relation: 



z=j (j{vBTH[5Br - (n^^r^E^n 

xexph^GL^r 5 ,^,^ 5 ^^)] 
= exp[-S cff (sl m) )]. (15) 

It is important to note that, when we adopt the Weyl symmetric formulation, 
we need to take into account the (Lorentz invariant) 5-function constraint 
among the dual gauge fields as 

x>i?r s = o. (i6) 
i=i 



As a result, the dual gauge symmetry of the DGL theory is kept in U(l) xU(l). 
In fact, although the DGL action (1) itself has the [U(l)] 3 dual gauge symme- 
try 

Xi i-> Xi exp(iai), B* eg i-> £f g - dai (i = 1, 2, 3), (17) 



this is finally reduced to U(l)xU(l) by the relation J2i=idcti = 0, which is 
implied by the constraint (16). Q 

Let us perform next the functional path integration over fields. For this aim, 
the quadratic term of the dual field strength in the action is evaluated as 



f (dBP + 2n * Cf } ) 2 = f {dB^f + 27r 2 /3 9 (j\ m \ A" 1 ^) , (18) 
where we use 



and the cross term vanishes. Next, inserting an identity for auxiliary fields £7j 
(1-form, i = 1,2,3), 

2 Here it may be interesting to comment that in this formulation we do not need an 
additional constraint for the phase of the monopole field such as Y2i=i Vi = const. 
In fact, such a requirement is already implicitly taken into account by Eq. (16). For 
instance, the special choice of the gauge, (unitary gauge) Qj = — rji (i = 1,2,3), 
leads to the corresponding constraint among the phase of the monopole field. 
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cxp 



1 /"p^exp - j^,-^) 



, (20) 



where, in general, Y%=i 7^ 0, an d rewriting the 5-function constraint for the 
dual gauge fields as a Fourier integral, the partition function takes the form 



z=jvk (i\vBTH\5B^-f Bl ^ (n^-^c 8 ) (n p ^) 

x exp [ - e( 2 ^« A " 1 ^ M ) + I (^ eS ) 2 + (#0 2 + 

• (21) 



i=l 



-i(Bf ng , ^) - i(sr g , ^ + dfc) - i(di^, ^) + _ „ 



2\2 



Here, A; is a 1-form auxiliary field. Notice that, due to (dr)™ 5 , Ej) = (r)™ 8 , <5-Ej), 
the integrations over the phase r/[ eg (i = 1, 2, 3) lead to the constraint SEi = 
which is resolved by introducing the Kalb-Ramond (KR) fields hi (2-form, 
i = 1,2, 3) as 



5[8Ei] = JvhiStfhi-fhiWEi-S*^], 



(22) 



where, in general, Y%=i hi ^ 0. Due to the hyper-gauge invariance h i— > h + dA 
with a 1-form field A, the hyper-gauge fixing functional also appears in order 
to avoid the overcounting in the integration over h. Then, we have 



Z = Jvk (j\VBTH[5Br - / B j) (iW 2 ) 
x (f[ VEiVhi5[5h t - f h i]5[Ei - 5 * hi]) 



\i=l 



x exp 



+\(E» j 2 Ei) - i(Bp*, Ei) - i{B?\ Ei + dk) + A(0 2 - v 2 f 



• (23) 



Now the integration over Ei (i — 1, 2, 3) leads to 
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z=fvk (jim es s[SBr s -fBi\j (um 2 ) (n^-/ ti 



x exp - jj^h (jt \ A- l jt } ) + | (rffiD 2 + {d^f 
+-(5 * hi, -U * hi) - t{Bf ns , 5* hi)- i(Bl es , 5*h t + dk) 



In order to evaluate the partition function further, we use the relations 



]T(ir s , s*hi)= 27T ^(a-m * s! m) , 5 * h 



(24) 



i=i 

3 



1 ( m ) > 



1=1 



i=l 



and 



(25) 



(5 * /ij, -g 5 * hi) = (dh u -pjdhi). 



By inserting an identity in form of a path-integration over j , 
const. = J Vf B i exp 



2£ 



(26) 



(27) 



and taking £ Bi = 1, we can integrate over £?[ eg as 

z = |m i^v^hiS^hi - f hi ]j 



x exp 



H — ^-(<5 * /ij + dk, A^(5 *hi + dk)) + (d<pi) 2 + ^-(dh h ^dhi) 



2Ki(hi, Sj m) ) - 2m(5K A^jT) + A(0 2 - v 2 ) 



1 -(m)\ 



(.2 „,2\2 



where 



(28) 



(5 * hi + dk, A-\5 * hi + dk)) = (h t ) 2 - (Shi, A~ l 5hi) + (kf 



(29) 
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Finally, the integration of k leads to a constant providing us with the expres- 
sion 



3 

nv^vhiSishi-u, 

\i=l 



x exp 



i=i I z Pg Z I J 9 



+ (dfaf + \{dK ^dhi) - 2m(h t , sj m) ) - 2m(6hi, A^f ) 



(30) 

As in our previous paper [39], we divide the action into three parts as 

S^SP + SW+SW], (31) 



where 



^-2n 2 (3 g (jt\A- 1 jt ) ), (32) 



SP = {dfaf + -(dhi, |^2-^| dht) + - v 2 ) 2 , (33) 

-27ri(/i i , £, M ) - 27ri(^, A-^ m) ). (34) 

Since these three actions have completely the same structures as in the SU(2) 
case [39] except for the labels of i, the further evaluation of these actions is 
straightforward. The first action provides the pure Coulombic potential for 
the static quark-antiquark system, since the color-electric currents interact via 
the Coulomb propagator. The second and third actions, and S^ 3 \ should 
be evaluated by taking into account the variation of the monopole field 0« 
in the vicinity of the center of the flux tube, which defines the "core" region 
of the flux tube. In fact, its classical profile shows the behavior fa = v at 
the distant region from the center, while it tends to vanish as fa — > on 
the Dirac string if ^ (if E-"^ = 0, the monopole field behaves as 

fa — v everywhere [29]). Effectively, the core region can be characterized by 
the inverse of the monopole mass as p < m^ 1 . It is interesting to note that 
the second action contributes only inside of the core since it vanishes in the 
region where fa = v , and its leading term contains the NG action with the 
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string tension <j corc [40]. At the same time, the third action in the core region 
contributes a boundary contribution Sf ore (j^). Then, the action for the core 
region is 



c 



= sr c ut ) )+v corc J d%4m, 



(35) 



where (a = 1, 2) parametrize the string world sheet Xi M (£«), and g(£i) is the 



determinant of the induced metric, g a b{^i) 



dxj M (&) dxj M (&) 



[40,41]. Here, the 



third action in the outside region p > m" 1 is responsible for the description 
of the "surface" of the flux tube. By inserting a suitable identity in form of a 
path- integral over the hyper-gauge fixing function fhi, 



const 



xcxp 



Av 2 



\ ([f h < + 4nvHDjt\ D-'A-'if, , + AmtiDj™]) 



,(36) 



where D = (A + m#) 1 is the propagator of the massive KR field, we can 
further integrate over the KR field as 



S, 



Ji 5 



D - A 



-i 



;( m ) 



(37) 



Now "| >m -i" means that a corresponding effective cutoff should be taken 
into account. One finds that the first term describes the interaction between 
the color-electric Dirac string via the propagator of the massive KR field. 
The derivative expansion of the KR propagator with respect to the covariant 
Laplacian on the string world sheet which is defined by 



9(Ci) 



(38) 



leads to the following form: 
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4vrV (Yt\DYt r 

= ^ surf / d%^~) 

s (m) 

+« surf / d% sfimg^i) (dj^i)) fat^m + O(Al), (39) 



where £^(£j) = ^J-jpj ^gg^ &X q^ is an antisymmetric tensor which deter- 
mines the orientation of the string world sheet xf™- 1 . The first term represents 
the NG action with the string tension 



t o m r + m 2 
a sml = 7cv 2 ln B \ x 



nv 2 In (l + /t 2 ) , 



(40) 



and the second term is the so-called rigidity term with the negative coefficient 



a 



surf 



71 V 



2 \m£ + m 2 B m 2 B/ 



7i(3 g k 2 

4 1 + K 2 



(< o), 



(41) 



where k = m x /mB is the GL parameter. 



Finally, by combining all contributions, the pure Coulomb term (32), the core 
term (35), and the surface term (37) with derivative expansion (39), we get 
the effective string action 



i=l 



^(jj m) ) + (a corc + a surf ) / d%^~) 

i 
i 

where the boundary (color-electric current) contribution is given by 

+2**p g (jlr\[D-A- i ]jir ) ) l - 

>m x 1 

Eqs. (42), (43) are the main results of this paper . 



, (42) 



(43) 



It has to be noted that the boundary term can be simplified only in the London 



limit m x — > oo. In this limit, the first term disappears, Sf oie (j 



Oh 



0, and 
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Fig. 1. The structure of the color-electric Dirac string in the color-magnetic repre- 
sentation for the mesonic (left) and the baryonic (right) quark systems. 

the third term becomes an "exact" Yukawa minus Coulomb term. Then, the 
cancellation of the Coulomb term takes place, and finally the boundary term 
is reduced into the ordinary Yukawa interaction as 



However, if the vacuum is not in the deep type-II limit, the boundary contribu- 
tion has, in principle, a complicated structure which can only be investigated 
numerically. 

Clearly, the explicit form of the effective string action depends on the quark 
system, that is to say, the combination of the color-electric charges. In the 
SU(3) case, one can consider not only the mesonic string, described by R-R, B- 
B and G-G systems, but also the baryonic string, as given by the combination 
of three color-electric charges R-B-G (see, Fig. 1). Remarkably, since the form 
of the action is now the Weyl symmetric one, which means invariant under the 
exchange of the color labels i, one can immediately conclude that all mesonic 
strings do not depend on the color. The baryonic system is constructed by 
connecting the three-colored mesonic strings, R-R, B-B, and G-G, at a certain 
junction, forming a Y-shaped structure. This is possible because the sum of 
all three color charges at the same point turns out to be zero. In Ref. [37], 
the static potential among three ends is calculated using the effective string 
action in the London limit and its short distance region is described by the 
Yukawa potential, which is consistent with our estimation in Eq. (44). The 
profile of the baryonic flux-tube solution in the DGL theory is investigated in 
Refs. [27,30], where also the role of the color-electric Dirac string inside the 
flux tube is clarified. 




(44) 
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5 Summary 



In this paper we have studied the effective string action of the color-electric 
flux tube in the dual Ginzburg-Landau (DGL) theory using the path-integral 
analysis in differential forms, thus extending our previous work for the U(l) 
dual Abelian Higgs (DAH) model [39] to the U(l)xU(l) case. By adopting 
the Weyl symmetric formulation, where the U(l)xU(l) symmetry is naturally 
extended to [U(l)] 3 , the treatment of the DGL theory is found to be rather 
simple, since the resulting form is the sum of three types of the U(l) DAH 
model. Taking advantage of this feature, we could similarly get the effective 
string action, which is expressed as the sum of the string action in the U(l) 
case. Then the interpretation of the action is also straightforward. The dif- 
ference is only the number of color-electric charges, which then allows us to 
describe not only the mesonic system but also the baryonic system. The latter 
is a typical system in the SU(3) gluodynamics. In the present paper, we do not 
describe the baryon systems explicitly, which will be investigated in a future 
analysis. 

Let us make some comments concerning the form of the effective string action 
in Eqs. (42) and (43). Analogously to the case of SU(2), the effect of the finite 
thickness of the flux tube leads us again to a modified Yukawa interaction 
as the boundary contribution of the open string, which turns over to the 
ordinary Yukawa one only in the London limit. The derivative expansion of 
the KR propagator with respect to the covariant Laplacian (38) leads to the 
NG action and the rigidity term as the leading and next-to-leading parts of 
the action. Such an expansion is clearly justified when the fluctuation of the 
string world sheet can be treated perturbatively. If not, one has to handle the 
expression of the first term in Eq. (37) directly. 

Finally, it is interesting to note that we can refer to the interaction among 
string world sheets of various colors without solving their dynamics. Since the 
sum of two types of color charges reduces to another type of antt-color charge, 
for instance as R + B — > G, one finds that if two strings composed of R-R and 
B-B interact with each other, this will be equivalent to a G-G string. Then, 
it is easy to imagine that the force between them must be attractive. 

One of the authors (Y.K.) is partially supported by the Ministry of Education, 
Science, Sports and Culture, Japan, Grant-in-Aid for Encouragement of Young 
Scientists (B), 14740161, 2002. We would like to thank D. Antonov and E.-M. 
Ilgenfritz for useful discussions. 



15 



References 

[I] G. 't Hooft, in High-Energy Physics. Proceedings of the EPS International 
Conference, Palermo, Italy, 23-28 June 1975, Vol. 2, edited by A. Zichichi, pp. 
1225-1249, Bologna, 1976. 

[2] S. Mandelstam, Phys. Rept. 23C, 245 (1976). 

[3] A. A. Abrikosov, Sov. Phys. JETP 5, 1174 (1957). 

[4] H. B. Nielsen and P. Olesen, Nucl. Phys. B61, 45 (1973). 

[5] Y. Nambu, Phys. Rev. D10, 4262 (1974). 

[6] A. S. Kronfeld, G. Schierholz, and U. J. Wiese, Nucl. Phys. B293, 461 (1987). 
[7] H. Shiba and T. Suzuki, Phys. Lett. B343, 315 (1995), |hep-lat/9406"oTo| . 
[8] H. Shiba and T. Suzuki, Phys. Lett. B351, 519 (1995), |iep-lat/9408004 . 
[9] Z. F. Ezawa and A. Iwazaki, Phys. Rev. D25, 2681 (1982). 
[10] Z. F. Ezawa and A. Iwazaki, Phys. Rev. D26, 631 (1982). 

[II] T. Suzuki and I. Yotsuyanagi, Phys. Rev. D42, 4257 (1990). 

[12] F. Brandstaeter, U. J. Wiese, and G. Schierholz, Phys. Lett. B272, 319 (1991). 
[13] O. Miyamura, Nucl. Phys. Proc. Suppl. 42, 538 (1995). 

[14] J. D. Stack, S. D. Neiman, and R. J. Wensley, Phys. Rev. D50, 3399 (1994), 
[hep-lat/9404014 . 

[15] O. Miyamura, Phys. Lett. B353, 91 (1995). 

[16] G. 't Hooft, Nucl. Phys. B190, 455 (1981). 

[17] S. Maedan and T. Suzuki, Prog. Theor. Phys. 81, 229 (1989). 

[18] T. Suzuki, Prog. Theor. Phys. 80, 929 (1988). 

[19] H. Suganum a, S. Sasaki, and H. Toki, Nucl. Phys. B435, 207 (1995), [hep] 
| ph/931235C . 

[20] S. Sasaki, H. Suganuma, and H. Toki, Prog. Theor. Phys. 94, 373 (1995). 

[21] V. Singh, D. A. Browne, and R. W. Haymaker, Phys. Lett. B306, 115 (1993), 
[hep-lat/9301004 . 

[22] Y.-c. Peng and R. W. Haymaker, Phys. Rev. D52, 3030 (1995), [hepj 
| lat/95030151 

[23] Y. Matsubara, S. Ejiri, and T. Suzuki, Nucl. Phys. Proc. Suppl. 34, 176 (1994), 
[hep-lat/9311061 . 



16 



[24] K. Schilling, G. S. Bali, and C. Schlichter, Nucl. Phys. Proc. Suppl. 73, 638 
(1999), |hep-lat/9809039| . 

[25] F. V. Gubarev, E.-M. Ilgenfritz, M. I. Polikarpov, and T. Suzuki, Phys. Lett. 
B468, 134 (1999), |hep-lat/9909099] 



[26] V. Bornyakov et al, Nucl. Phys. Proc. Suppl. 106, 634 (2002), |h^pdat/0111042 



[27] S. Kamizawa, Y. Matsubara, H. Shiba, and T. Suzuki, Nucl. Phys. B389, 563 
(1993). 



[28] Y. Koma, H. Suganuma, and H. Toki, Phys. Rev. D60, 074024 (1999), [hepj 



ph/9902441. 



[29] Y. Koma and H. Toki, Phys. Rev. D62, 054027 (2000), |hep-ph/0004T77 . 



[30] Y. Koma, E.-M. Ilgenfritz, T. Suzuki, and H. Toki, Phys. Rev. D64, 014015 
(2001), |hep-ph/0011165| . 



[31] P. Orland, Nucl. Phys. B428, 221 (1994), jhep-th/9404l"4C| . 



[32] M. Sato and S. Yahikozawa, Nucl. Phys. B436, 100 (1995), |hep-th/940620^ . 

[33] E. T. Akhmedov, M. N. Chernodub, M. I. Polikarpov, and M. A. Zubkov, Phys. 
Rev. D53, 2087 (1996), |hep-th/9505070 . 



[34] D. Antonov and D. Ebert, Eur. Phys. J. C8, 343 (1999), [hep-th/9806153 . 
[35] A. M. Polyakov, Nucl. Phys. B268, 406 (1986). 
[36] H. Kleinert, Phys. Lett. B174, 335 (1986). 



[37] M. N. Chernodub and D. A. Komarov, JETP Lett. 68, 117 (1998), jaej* 
th/9809l53| . 

[38] D. Antonov and D. Ebert, Phys. Lett. B444, 208 (1998), |hep-th/9809018j 

[39] Y. Koma, M. Koma, D. Ebert, and H. Toki, JHEP 08, 047 (2002), [hep 
| th/010813^ . 

[40] D. Foerster, Nucl. Phys. B81, 84 (1974). 

[41] R. L. Davis and E. P. S. Shellard, Phys. Lett. B214, 219 (1988). 



17 



